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commodity cultivation of plants is carried out. During formulation of 
the model, inter alia, models described in literature and results of 
experimental research were used. A developed mathematical model 
was implemented in MATLAB/Simulink software and simulations 


aa carried out with a computer model were used for carrying out 
microclimate graphical and statistical validation of a model. Analysis of simulation 
mathematical model results allows making statement of logical correctness of the 
heat and mass exchange developed model and makes it possible to determine critical points of 


failure to adjust the model. A simplification degree of the developed 
heat exchange model influences its precision. In order to use the 
developed model e.g. for the control purposes, it requires to be more 
detailed. 


Introduction 


The interest in the greenhouse climate modeling dates back to the commodity crops 
under cover. Greenhouse climate models are the main tool helping to control thermal and 
humidity parameters inside the greenhouse. Knowledge of them supports decision-making 
related to cultivation, and also allows the use of modern, complex microclimate control 
algorithms. In the literature we can find many greenhouse climate models, they are both 
static and dynamic mathematical models, as well as „black box” models based on artificial 
intelligence methods. A review of these models can be found in the works (Boaventura 
Cunha, 2003; Raczek, 2012). Among the most important in this respect are the works of 
Bot (1983) and Jong (1990). Most greenhouse climate models are formulated for the 
experimental facilities, equipped with many additional sensors. Often these models are 
deliberately simplified. Some processes are omitted in the description in order to investigate 
a phenomenon which is interesting for the researchers (e.g. studies on ventilation are 
carried out in the facility without plants). The greenhouse climate models are formulated 
and validated taking into account: the crop type and development phase of crops, region 
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and weather conditions; the structure and type of a greenhouse, and the operation of 
ventilation equipment as well. Therefore, it is not easy to directly extrapolate these models 
to differently built greenhouses located elsewhere. 

The object of the research in this paper was modern Venlo type greenhouse located in 
Różanki in Lubuskie Voivodeshop, where on 5.9 hectares tomato cultivation was carried 
out. For the greenhouse a mathematical model of heat and mass transfer processes in the 
indoor air was developed which is presented below. In March and April 2011 temperature 
and relative humidity of air inside the greenhouse was recorded by a climate computer 
every 15 minutes, which in the model constitute the output signals, and control signals 
(input in the model): in the form of the heating pipe temperature and the degree of opening 
vents. Environmental parameters of the greenhouse constituted interference: temperature 
and humidity of the outside air, radiation and wind speed. The weather data were recorded 
by a weather station. 

The objective of the paper was to draw up a dynamic, mathematical model of heat and 
mass exchange in air inside a big-size greenhouse, where a commodity cultivation of plants 
is carried out when we are provided with data on the climate parameters normally collected 
by a climate computer and a weather station. 

The scope of work related to the greenhouse climate modeling included: 

— adaptation of literature models to the test object, by taking into account characteristic 
dimensions of the tested greenhouse and its technical equipment, 

— draw up a computer model and perform simulations, 

— graphical and statistical validation of the resulting model of the process of heat 
exchange and mass transfer in the air inside the greenhouse using the results of 
experimental studies. 


Mathematical model of greenhouse microclimate 


Microclimate in the greenhouse is the result of combination of complex mechanisms 
involving processes of heat and mass transfer occurring in the greenhouse, and the 
processes between the interior of the greenhouse and the surroundings. Processes occurring 
in the greenhouse are highly non-linear, related to each other. During the formulation of a 
mathematical model of these processes, equations of heat and mass balances for the air 
inside the greenhouse should be formulated. (Wachowicz, 2006). 

To develop a model the following simplifying assumptions were made: 

— a greenhouse is treated as a perfectly mixed tank, i.e. the analyzed air parameters have 
the same value in the entire volume of the greenhouse, 

— due to the lack of empirical data no effect of energy screens was taken into account, 

— impact of the ground on the heat and mass transfer was neglected, because in the tested 
facility, cultivation was carried out on coconut fiber mats covered with white foil, 

— evaporation from the greenhouse cover and crops was neglected, because in modern 
greenhouses a condensate is drained from the cover, and when cultivation is maintained 
properly, retting on plants should not take place. 
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Heat transfer model 


Considering the equipment of the object, changes in air temperature inside the 
greenhouse can be presented in the form of the heat balance equation: 


Wa = = mA (Oriona + O soraw F O mdiagja = Owent m O ransp + Orona ) (K ` s”) (1) 
where: 

Trew — air temperature inside the greenhouse, (K) 

Pwew — air density inside the greenhouse, (kg:m°) 

Önen — specific heat of air inside the greenhouse, (Jkg"-K") 

Vs. — greenhouse volume, (m*) 

Qostona  — heat flux exchanged between the interior and the cover, (J ‘s') 


Q;.grzewcz — heat flux from heating pipes, (J s") 
Qradiacja — heat flux supplied from solar radiation, (J s") 


Qwent — heat flux exchanged through the ventilation, (J s") 
Oransp — heat flux exchanged by plant transpiration, (J ‘s') 
Okond — heat flux supplied by condensation, (J-s"') 


Heat balance components are heat fluxes supplied to the greenhouse: due to 
condensation of water vapor on the cover, from heating pipes, from solar radiation and heat 
fluxes exchanged during ventilation, and as a result of plant transpiration. Heat transfer 
through the cover may take place in both directions. 

Most of convective heat fluxes exchanging between different parts of the greenhouse 
and the air inside the greenhouse depend on the heat transfer coefficients and the 
temperature difference between the elements surface and the air. These fluxes describe the 
following equations: 

— for convective heat exchange between the air inside the greenhouse and the cover: 


Qoslona = Qoslona A slona (Rów = Trew ) (J-s') (2) 


where: 
Oostona — heat transfer coefficient through the cover, (W-m~-K”) 
Aoslona — cover surface, (m°) 
Tostona — cover temperature, (K) 

— for convective heat exchange between the air and the heating pipes: 


-1 
0. grzewczy — as. grzewczy A; grzewczy (7 s.grzewczy | T wew ) (J 'S ) (3) 


where: 
Gs. grzewczy — heat transfer coefficient of heating pipes, (W.m”-K"') 
As erzewezy — heating pipes surface, (mô 
T, grzewczy — heating pipes temperature, (K) 
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— for convective heat exchange through ventilation: 


Owent = Avent Awent (Tew m dY (J 8") (4) 
where: 
Owen — heat transfer coefficient between the air inside the greenhouse and its 
environment, (W-m”-K”') 
Avent — vent sufrace, (m”) 
Tzw — Outside air temperature, (K) 


Heat flux supplied from solar radiation is expressed by following simplified relation 
(Tap, 2000): 


Q ,adiacja = A,,0,7Rad (Js) (5) 


where: 
Rad —solar radiation density, (W-m”) 


From analysis of the results of simulation studies of the greenhouse follows that, in this 
case value 0.7 in the equation (5) is too high and was lowered to the level of 0.55. 

Determination of the heat transfer coefficients a is a complex task. The reason for this is 
primarily a large number of time-varying factors that affect the value of coefficients. In this 
study, the heat transfer coefficients are calculated on the basis reported in the literature 
information models (Bot, 1983; Zwart, 1996), i.e. experimentally determined relationship 
between the temperature difference of the elements and the air temperature in the 
greenhouse and characteristic values. These coefficients are expressed by the following 
formulas: 


— heat transfer coefficient between the air and the inner side of the cover (Zwart, 1996): 


Qosiona 7 1,7(cos oP” ` (T =, |= (W-m*.K"') (6) 


wew oslona 


where: 
p — Slope of the roof, (°) 
— heat transfer coefficient between the air and the upper heating pipes, which are located 
above the plants (Bot, 1983): 
Q s grzewczy | 7 1,28 a... zm = yet p” (W-m*.K"') (7) 


s.grzewczy 1" 


— heat transfer coefficient between the air and the lower and vegetative heating pipes, 
when the pipes have a diameter of 51 mm and are located in the area of plants growth 
and under the cultivation gutters (Bot, 1983): 
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32 -2 œ -1 
Q s grzewczy 2 = 1,99 A, grzewczy 2° (r s.grzewczy 2 T wew p (W-m K ) (8) 


— heat transfer coefficient between the air inside the greenhouse and its surroundings 
through natural ventilation (Kurpaska, 2007): 
Q went = P wew i C wew ` D (W-m*.K") (9) 


went 


where: 
Pwent — air flow rate through ventilators, (m-s) 


taking into account density and specific heat of the air inside the greenhouse, and air flow 
rate through ventilators. The simplified relation presented in the work by Tap (2000) is used 
in the present work to calculate the air flow rate through ventilators 


D vent (Ee K, | (m-s") (10) 
1+4:K, 
where 
K, — opening degree of ventilators on the windward side, (%) 
K, — opening degree of ventilators on the leeward side, (%) 
Way — outside wind speed, (m-s”) 


o,y,l,EĘ,V — constants, the value of which amounts to: 6=7.1708-10°(%"'), y=0.0156(%"), 
C=2.7060-10°(-), €=6.3233-10°(%"'), y=7.4-10° (m:s"') 


Air flow rate depends on the degree of opening vents on the windward side K, and lee 
side K, as well as on the outside wind speed vew. 

The last two heat fluxes Qrransy and Okona included in the heat balance are latent heat. 
One of them is heat O, lost through transpiration of plants (11), when due to plant 
metabolism the water passes from the liquid phase into the gaseous one. 


Wag = Fo Aliscie M transp (J-kg”') (l 1) 
where: 
ro — heat of vaporization or condensation, (J-kg") 
Anicie  — leaf area, (m°) 


Miransp — Mass flux absorbed by plants through transpiration referred to m” of the 
greenhouse area (KgpąryM +8") 


The second source of the latent heat is the condensation on the cover. By condensation 
of water vapor, heat is released to the environment: 


Okona = ToÁoslonaM kond (J-kg”) (12) 


where: 
. Fl 
Mbkond — mass flux from the condensation (kgpary.M *S ) 
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Mass exchange model 


Water vapor content of the air inside the greenhouse is an important environmental 
parameter, which determines the proper development of a crop. It is also used to assess the 
risk of diseases and undesirable pest development. 

A model describing the changes in the water vapor content of the air inside the 
greenhouse is based on the mass balance equation. The primary source of water vapor in 
the balance equation is plant transpiration. Water vapor in the air inside the greenhouse is 
reduced as a result of condensation on the inner side of the cover. Mass transfer due to 
ventilation can take place in both directions depending on the conditions inside and outside 
the greenhouse. Mass balance equation takes the following form: 


Vee A wew 
"ak a > M ransp M kond M went (13) 
SZ 
where: 
As- — greenhouse area, (m°) 
Jfwew — Water content of the air inside the greenhouse, (gm ) 


a: 2-1 
Mwen — mass flux due to ventilation, (kgpay-m™~-s” ) 


In this study to determine the transpiration of the plants a regression model presented in 
the work (Kurpaska, 2006) was applied. It was developed to determine the water demand of 
greenhouse tomatoes, taking into account controllable factors of the surrounding climate, 
i.e. solar radiation Rad, air temperature inside the greenhouse 7,,„ and vapor pressure 
deficiency VPD. 


M ransp = — (0.0025 -Rad + 0,098 -T,,,,, — 0,143 -VPD + 0,05) (g-m™-s) (14) 
where: 
VPD -vapor pressure deficit, (Pa) 


In order to simplify the developed mass exchange model it was assumed that the 
temperature of inner side of the cover is calculated according to the following equation 
(Kurpaska, 2007): 


dostana = 0,47, 


zew 


+0,6T,,„ (K) (15) 


where: 
Tostona — cover temperature, (K) 
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Whereas condensation of water vapor on the inner side of the cover is expressed by 
the relation (Tap, 2000): 


my 
M rona =m; Trew a L tna | z Tasan gdy ba 7 F vax MAF (g-m”-s) (16) 
M kond T 0 gdy Ja < La wew 
where: 
m. . 
My IT, wee ae * — mass transfer coefficient, 
m;im,  — fixed parameters of the mass transfer coefficient, 
Jmaxwew  — the maximum water content of the air inside the greenhouse (gm”). 


Reducing or increasing (to a lesser degree) the water vapor content of the air inside the 
greenhouse is also carried out through ventilation. The amount of water vapor removed 
from the greenhouse is described by the following formula: 


ae = ® seni (Few > lx) (gm”*.s) (17) 


where: 
few  — Water content in the outside air, (m°) 


The air flow rate through ventilators is calculated from the relationship (10) presented in 
the part of the work on heat balance. 


Computer model and simulation results 


Models of heat and mass transfer in the air inside the greenhouse, which were described 
earlier, were implemented in MATLAB/Simulink software. Figure 1 shows a schematic 
representation of the developed computer model of heat and mass transfer in the air inside 
the greenhouse. 

The developed computer model was used to perform simulations, which led to the 
conduct of daily changes in temperature and humidity inside the greenhouse. Waveforms of 
input variables of the simulation model shown in Figure 1 were determined by experimental 
investigations. Figure 2 and figure 3 show exemplary results of these simulations for four 
consecutive days, as well as the values of the analyzed parameters obtained in the course of 
experimental tests in the greenhouse. The simulations carried out with a computer model 
were used for carrying out graphical and statistical validation of the model. 
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Heat exchange model 


transpiration 
ventilation 


ventilation 


Mass exchange model 


Figure 1. Schematic representation of the computer model of heat and mass exchange in air 
inside greenhouse. Symbols: Wywew— relative humidity of air inside greenhouse (%), Wzew— 
relative humidity of air inside greenhouse (%), Tsym — air temperature inside 
greenhouse obtained as a result of computer simulations (°C), W sym — humidity of air inside 


> $ > =3 
greenhouse obtained as a result of computer simulations (kZpary M” powietrz 


experim. data 
simulations 


"0 1000 2000 3000 4000 5000 
Time: 1440 = 24h 


Figure 2. Air temperature changes inside greenhouse obtained during measurements and 
simulation research 
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Analysis of simulation results allows statement of logical correctness of the developed 
model. It also enables the identification of critical points of failure to adjust the model. It 
was found out that in this case the moment is the sunset, marked on the charts with 
a vertical line 9= =®. This indicates the strong correlation and sensitivity of microclimate 
inside the greenhouse on the amount of heat supplied to the object due to sun radiation. 
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Figure 3. Changes of air humidity inside greenhouse obtained as a result of measurements 
and computer simulations 


For statistical verification of the developed model a determination rate was used (Makać 
and Urbanek, 2010): 


re Al (18) 


where: 
Yi — are the values of the selected features obtained from measurements, 
Y; — the corresponding values determined from the model. 


If the value of determination rate is closer to 1, the better is consistency of data from the 
model with empirical data. For the heat exchange model the determination rate (for 
example four days presented in fig. 2 and 3) was 0.87. Whereas for mass exchange model 
the determination rate was 0.97. 
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Conclusions 


The accuracy of the heat transfer model is mainly affected by the simplification degree 
and by omitting other heat fluxes in the model. This applies in particular to the heat 
exchanged between solid elements of the greenhouse and the air inside the greenhouse by 
radiation. None of these heat fluxes in the developed model may explain the discrepancy 
between the data obtained during the simulation and data obtained from experimental 
investigations in the period just after sunset. Whereas a better fit of the mass transfer model 
may result, inter alia, from smaller number of process variables, compared to the amount of 
variables involved in the heat exchange process. The developed mathematical model and 
a computer model created on this basis provide valuable information about the complex 
system, which is a greenhouse. The ability to perform a variety of simulation enables, 
among others, to analyze the sensitivity of the system to individual parameters changes. In 
order to use the developed model e.g. for control purposes, it requires to be more detailed. 
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MODEL PROCESU WYMIANY CIEPŁA I MASY 
W POWIETRZU WEWNĄTRZ SZKLARNI 


Streszczenie. Celem pracy było opracowanie matematycznego modelu wymiany ciepła 1 masy 
w powietrzu wewnątrz wielkogabarytowej szklarni, w której prowadzona jest towarowa uprawa 
roślin. Podczas formułowania modelu wykorzystano m.in. modele opisane w literaturze i wyniki 
badań eksperymentalnych. Opracowany model matematyczny został zaimplementowany do programu 
MATLAB/Simulink, a symulacje przeprowadzone z udziałem modelu komputerowego wykorzy- 
stano do przeprowadzenia graficznej i statystycznej walidacji modelu. Analiza wyników symulacji 
pozwala na stwierdzenie logicznej poprawności opracowanego modelu, a także umożliwia określenie 
punktów krytycznych niedopasowania modelu. Na dokładność opracowanego modelu wymiany 
ciepła wpływa przede wszystkim stopień jego uproszczenia. Aby opracowany model mógł być 
wykorzystany, np. do celów sterowniczych, wymaga większego uszczegółowienia. 


Słowa kluczowe: szklarnia, mikroklimat, model matematyczny, wymiana ciepła i masy 
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